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Abstract 

The critical behavior of a complex iV-component order parameter Ginzburg-Landau 
model with isotropic and cubic interactions describing antiferromagnetic and structural 
phase transitions in certain crystals with complicated ordering is studied in the frame- 
work of the four-loop renormalization group (RG) approach in (4 — e) dimensions. By 
using dimensional regularization and the minimal subtraction scheme, the perturbative 
expansions for RG functions are deduced and resummed by the Borel-Leroy transfor- 
mation combined with a conformal mapping. Investigation of the global structure of 
RG flows for the physically significant cases N = 2 and N = 3 shows that the model 
has an anisotropic stable fixed point governing the continuous phase transitions with 
new critical exponents. This is supported by the estimate of the critical dimensionality 
N c = 1.445(20) obtained from six loops via the exact relation N c = ^n c established for 
the complex and real hypercubic models. 
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I Introduction 



In this paper we investigate the critical properties of phase transitions in certain antifer- 
romagnets involving an increase of the unit cell in one or more directions at the critical 
temperature. They are known to be described by a generalized 2n-component (n > 2) 
Ginzburg-Landau model with three independent quartic terms 

„ i 2n 2n „ 

8=1 ' 1=1 

4 - i=l 4 - i=l 

associated with isotropic, cubic, and tetragonal interactions, respectively JTJ. Here <pi is 
the real vector order parameter in D = 4 — e dimensions and is proportional to the 
deviation from the mean-field transition point. According to the universality hypothesis, 
the critical phenomena in various physical systems must not depend on microscopic details, 
but are determined solely by the spatial dimensionality of the system, the interaction range 
as well as the symmetry and the dimensionality of the order parameter field. The anomalous 
behavior of the thermodynamic quantities and the correlation radius r c are expressed by 
power laws in the critical region. 

When n — 2, Hamiltonian (|l|) describes the antiferromagnetic phase transitions in TbAu2 
and DyC2 and the structural phase transition in the NbG"2 crystal. The phase transitions 
in the helical magnets Tb, Dy, and Ho belong to the same class of universality [[|. Another 
physically important case n = 3 is relevant to the antiferromagnetic phase transitions in 
such substances as K^IrCle, TbD2, and Nd. The phase transition in antiferromagnet MnS2 
belongs to the same class of universality 0. All these phase transitions are predicted to be 
of second orderQ |2|, E|. This is confirmed by several experiments (see Ref. || and references 
therein). However the experimental data were insufficiently accurate to provide reliable val- 
ues of critical exponents and the obtained estimates |7|, ^ |S| were found to differ significantly 
from the theoretically expected numbers. Such a distinction should be interpreted as a man- 
ifestation of influence of various possible effects in the real substances such as the presence 
of defects. 

For the first time, the magnetic and structural phase transitions described by model 
([[]) were studied in the framework of the renormalization group (RG) by Mukamel and 
Krinsky within the lowest orders in e [Q. A three-dimensionally stable fixed point (FP) with 
coordinates u* > 0, v* = z* > was predicted^. That point was shown to determine a 
new universality class with a specific set of critical exponents. However, for the physically 
important case n — 2, the critical exponents of the unique point turned out to be exactly 
the same as those of the 0(4)-symmetric one. 

3 An interesting type of multisublattice antiferromagnets, such as MnO, CoO, FeO, and NiO, was studied 
in Ref. || ||, It was shown, in the leading orders in e, that the phase transitions in these substances are 
of first order. 

4 Following Mukamel we call this point "unique". 
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For the years, an alternative analysis of critical behavior of the model, the RG approach 
in three dimensions, was carried out within the two- and three-loop approximations (TO, 



11]. Those investigations gave the same qualitative predictions: the unique stable FP does 
exist on the 3D RG flow diagram. By using different resummation procedures, the critical 
exponents computed at this point proved to be close to the exponents of the Bose FP 
(it = 0, v = z > 0) rather than the isotropic (O(n)-symmetric; u > 0, v = z = 0) one. It 
was also shown that the unique and Bose FPs are very close to each other, so that they may 



interchange their stability in the next orders of RG approximation [11 



Recently, the critical properties of the model were analyzed in third order in e [ p"2| . [TjJ . 
Investigation of the FP stability and calculation of the critical dimensionality n c of the 
order parameter, separating two different regimes of critical behavior^, confirmed that model 
(0) has the unique stable FP at n = 2 and n — 3. However, the twofold degeneracy of 
the stability matrix eigenvalues at the one-loop level was observed for this FP [13]. That 



degeneracy was shown to cause a substantial decrease of the accuracy expected within the 
three-loop approximation and powers of yfe to appear in the expansions^. So, computational 
difficulties were shown to grow faster than the amount of essential information one may 
extract from high-loop approximations. That resulted in the conclusion that the e-expansion 
method is not quite effective for the given model. 

Another problem associated with model (H) is the question whether the unique FP is 
really stable in 3D, thus leading to a new class of universal behavior, or its stability is only 
an effect of insufficient accuracy of the RG approximations used. Indeed, there are general 
nonperturbative theoretical arguments indicating that the only stable FP in 3D may be the 



Bose one and the phase transitions of interest should be governed by that stable FP [15 



The point is that when v — z, the model (]]]) describes n interacting Bose systems. As was 



shown by Sak fL6fl , the interaction term can be represented as the product of the energy 
operators of various two-component subsystems. It was also found that one (the smallest) of 
the eigenvalue exponents characterizing the evolution of this term under the renormalization 
group in a neighborhood of the Bose FP is proportional to the specific heat exponent a. Since 
a is believed to be negative at that point, and that is supported experimentally [O and 



theoretically [|T^, [19], |2T]], the interaction is irrelevant. Consequently, the Bose FP should 
be stable in three dimensions. However, up to now, this conclusion found no confirmation 
within the RG approach, so the model (0) poses a certain challenge to the perturbative 
methods in the theory of critical phenomena. It is therefore highly desirable to extend 
already known e expansions for the stability matrix eigenvalues, critical exponents and the 
critical dimensionality in order to apply more sophisticated resummation techniques to longer 
expansions. 

Thus, in the present work we firstly avoid the problem of the eigenvalues degeneracy in 
model (|TJ) by analyzing the critical behavior of the equivalent complex iV-component order 

5 Whcn n > n c the unique FP is stable in 3D while for n < n c the stable FP is the isotropic one. 
6 Similar phenomenon was observed earlier in studying the impure Ising model (see Ref. 14 ). Half-integer 
powers in e arising in that model have different origin but also lead to the loss of accuracy. 
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parameter model with the effective Hamiltonian 



r 1 N u N v N 

H = J d ° x [2 EKM + vViW*) + + E itoMW 

j=l ' i,j=l ' i=l 



(2) 



comprising the isotropic and cubic interactions^. Note that this Hamiltonian also determines 



the real hypercubic model |^|] relevant to the magnetic and structural phase transitions in 
a cubic crystal if ijji is thought to be the real n-component order parameter. The model (0) 
comes out exactly from model ([!]) at Vq = z and it is free from the eigenvalues degeneracy. 
Secondly, we examine the existence of an anisotropic stable FP in model @ on the basis 
of the higher-order e expansion. Namely, using dimensional regularization and the minimal 
subtraction scheme ||24|| , we derive the four- loop RG functions as power series in e and 
analyze the FP stability. For the first time, we give numerical estimates for the stability 
matrix eigenvalues from e expansions using the Borel-Leroy transformation with a conformal 
This allows us to carry out the careful analysis of the stability of all the 

\n c between the critical (marginal) 



25, 26 



mapping 

FPs of the model. We establish the exact relation N r 



spin dimensionalities of the complex and the real hypercubic models and obtain the estimate 
N c = 1.445(20) using six- loop results of Ref. ||27|| . We show that the anisotropic (complex 
cubic; u ^ 0, v ^ 0) stable FP of model (|2|), being a counterpart of the unique point in 
model flip, does exist on 3D RG flow diagram at N > N c . For this stable FP we give more 



accurate critical exponents estimates in comparison with the previous three-loop results [13 
by applying the summation technique of Ref. |28[ to the longer series. 



II Four-loop e-expansions, resummation and the fixed 
point stability 

To deduce RG expansions for the /3-functions and critical exponents one needs to calculate a 
set of Feynman graphs, each of them comprising three factors: the combinatorial coefficient, 
the result of tensor contractions and the integral value associated with the diagram. The 
combinatorial factors and the values of integrals were found earlier in Ref. p6j . To evaluate 
the tensor contractions for the vertex and mass diagrams, the tensors G\ and G2 correspond- 
ing to the isotropic and cubic interactions in Hamiltonian (0) are introduced; they have a 
simple symmetrized form: 

2ijki = Oijdikdii , 

where {a, j3, fx, u} — 1,2 and {i, j, k, 1} — 1, . . . , TV. 

7 The model with the complex vector order parameter was considered by Dzyaloshinskii in studying 
the phase transitions in DyC 2 , TbAu 2 (N = 2) and TbD 2 , MnS 2 , and Nd (N = 3). 
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Further, normalizing conditions on the renormalized one-particle irreducible inverse Green 
functions 1^ and vertices given by corresponding Feynman diagrams must be imposed. 
Within the massless theory they are as follows 





1 



\i e u 

jJL e V 



(3) 



with one more condition on the \ip\ 2 insertion 

,(1,2), 



1 9 



(4) 



Here m, u, and v are the renormalized mass and dimensionless coupling constants, fi is an 
arbitrary mass parameter introduced for dimensional regularization. The vertices and 
ri 4 ) are connected with the vertex function without external lines 



p(4) a0uv 
ijkl 



p(4) (~i a(3fiu p(4) (~i a/3 fiu 



^From renormalization conditions ([|) and the expansions for the renormalization con- 
stants Z^p, Z u , Z v , and Z\m2 may be obtained. These constants relate the bare mass m and 
coupling constants Uo, Vq of the Hamiltonian (0) to the corresponding physical parameters: 



Z,i, 



m 



7 2 

Zj m m 



u 



e 



"0 



(5) 



So, with relations ([5]) taken into account, the /^-functions and critical exponents can be 
calculated via the formulas 



d In up n | d In up a 
du ^ u ~r dv & 



d In vp o I d In tip a 
du ^ dv 



T](U,V) 



= -e , 
= -e , 

a _|_ ydln_2^p_n 
du r J u ~r ^ g v /J V , 



(6) 



(7) 



9m ^ m 9u 

where /3 9 = grf^j , # = {^,^}- The critical exponents r/ and t? 2 are found by substituting 
zeros of the /3-functions into expressions (0). The critical exponent 7 is calculated through 
the known scaling relation 7 -1 = 1 + 



The four-loop e expansions for the /3-functions of model 
Ref. ra. They read 



eu — u — 



N + 4 



1 



(TV + 4) 2 
4 



were recently obtained in 
3w 3 (3iV + 7) + 44u 2 v + llW 1 



~ (AT + 4)3 [x( 48C(3)(5iV + U) + 33Ar2 + 461iV + 74 °) 

/ \ U 2 V 2 / \ 

+ m 3 v(384C(3) + 79iV + 659) + — (288C(3) + 3iV + 1078) + 14 W 

~ (N I 4)4 [^(" 48 C(3)(63iV 2 + 382iV + 583) + 144C(4)(5iV 2 
+ 31iV + 44) - 480C(5)(4A^ 2 + 55N + 93) + 5iV 3 - 3160iV 2 
- 20114iV - 24581) - ^-^(l2C(3)(3iV 2 + 2767V + 1214) - 36C(4) 



x 



3 

(197V + 85) + C(5)(2400/V + 23040) - 28iV 2 + 3957/V + 15967) 



u 3 v 2 



(72C(3)(19iV + 426) - 4032C(4) + 39840C(5) + 1302iV + 46447) 



+ — ^— (60C(3) (N - 84) - 792C(4) - 4800C(5) - 125/V - 12809) 

4 

HI) / \ -i 

_(400C(3) + 768C(4) + 385l)] , 

V = ev- -J—Oouv + 5v 2 ) + — 1 — \u 2 v(5N + 41) + 80uv 2 + 30v' 
7V + 4 V ' (N + 4) 2 L v ' 

^4p [~2~ ( 96 C( 3 ) + 7) - 13iV 2 + 184iV + 821 



+ 



(N + 

(4032C(3) + 59iV + 5183) + w 3 (768C(3) + 1093) + V — (384C(3) + 617 



u 2 v 2 / ,„ s r \ ,/ \ v 4 



- (Ar | 4)4 [^(48C(3)(iV 3 - 12iV 2 - 140iV - 567) + 144C(4)(2iV 2 
+ 17iV + 45) - 3360C(5)(3iV + 13) - 29iV 3 - 28iV 2 - 6958iV - 19679) 



+ ^(l2C(3)(9iV 2 - 5917V - 7028) + C(4)(3528iV + 21240) - 480C(5) 

x (107V + 287) + 61iV 2 - 5173/V - 66764) - ^(l800C(3)(/V + 62) 

- 144C(4)(8iV + 203) + 172800C(5) + 56iV + 9370l) - 4wt> 4 (5090C(3) 

- 1296C(4) + 7600C(5) + 4503) + y (-8224((3) + 1920C(4) 

- 12160C(5) - 7975)" 

where C(3), C(4), and £(5) are the Riemann ( functions. For the critical exponents we fi 

u 2 {N+ 1) + 4uv + 2v 2 
-(N + 1)(N + 4) + 3u 2 v(N + 4) + 15uv 2 + 5v d 



1 



(N + 4) 2 

1 r v* 
(iV + 4) 3 



+ 



(iV + 4) 4 



■-u*(N + 1) (iV 2 -9N- 25) - Wu 3 v(n 2 

165 



- 9N - 25 ) + 495mV + 330uv a + 



(10) 



n-2 



1 



7V + 4 
1 

(JV + 4) 3 



2u(iV + 1) +Av 
3 



(iV + 4) 5 



6w 2 (iV + l) + 24m; + 12t/ 



u 



31N 2 + UQN + 115 + 3u 2 v(31N + 115) 



+ 3uv 2 (N + 145) + 146v E 



+ 



(N + 4) 4 



(l2((3)(3N 3 + 8N 2 



+ 22N + 17) + 36C(4)(5iV 2 + 167V + 11) - 4iV 3 + 977iV 2 

+ 30167V + 2035) + -^-(l2C(3)(3iV 2 + 5iV + 17) + C(4)(180iV 

+ 396) - 4iV 2 + 981iV + 2035) + Au 2 v 2 (((3) (18 JV + 132) 

2uv 3 



+ C(4)(18JV + 270) + 1157V + 1391 



60C(3)(iV- li; 



1152C(4) - 1257V - 5899) + 4w 4 (25C(3) + 48((4) + 251 



(11) 
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1 



(TV + 4) 
1 



u(N+ l) + 2v 



u 2 {N + 1) + Auv + 2v 2 



(N + 4) 2 

3u 3 (N + 1)(UN + 39) + 18w 2 t;( lliV + 39) 
1 



: (24C(3)(iV+ 1)(3/V 2 



4(iV + 4) 3 

+ 10m; 2 (jV + 89) + 300t; 3 ] + — - - u 

+ 5iV + 17) + 72((4)(iV + l)(5iV + 11) - 5(N + 1)(N 2 - 399N 

- 820)) + 8w\(24C(3)(3iV 2 + 5N + 17) + C(4)(360iV + 792) 

- 5(iV 2 - 399iV - 820)) + 6uV(48C(3)(3JV + 22) 

+ 144C(4)(A^+ 15) + 953A^+ 11227)-2^ 3 (240C(3)(iV- 11) 

- 4608C(4) - 5157V - 23845) + 12w 4 (l00C(3) + 192C(4) + 1015)] . (12) 

^From the system of equations /3 u (u*,v*) = 0, v (u*,v*) = 0, one can calculate all the 
FPs of model (|2|), which are given by the power series in e: 



u* = u*(e) = u k£ k > 

k=l 



V = V 



k=l 



There exist four FPs, one of them (Gaussian) is trivial: 
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1. Gaussian FP 

u* = v* = . 

2. Isotropic or 0(iV)-symmetric FP 

u * = e+^(3^ + 7)e 2 - i( ^(48C(3)(iV + 4)(5iV+ll) 

+ 33 N 3 - 55iV 2 - 440iV - 568) e 3 

+ i2(w( 48 C(3)(^V + 4)(63iV 3 - 417V 2 - 9497V - 1133) 

- 144C(4)(iV + 4) 3 (5A r + 11) +480C(5)(iV + 4) 2 (4iV 2 
+ 55iV + 93) - 5iV 5 - 1335iV 4 - 2(698iV 3 - 3299iV 2 

- 9666iV-8278)e 4 , 
v* = . 

3. Bose FP 

u* = , 

v* = |(iV + 4)£+|(iV + 4)e 2 + T ^(iV + 4)(l03 - 384C(3))5 3 
+ 62W ( N + 4 ) (~3296C(3) - 1920C(4) + 12160C(5) + 265)e 4 . 

4. Complex cubic FP 

u* = ^- 4 (N + A)e+ T ^ w (N + A)(70N 2 -205N + 139)e" 

+ l(5lhw( N + 4 ) (48C(3)(5iV - 4) (64iV 3 - 188iV 2 + 1517V 

- 23) - 6370iV 4 - 24149iV 3 + 1447197V 2 - 197208iV + 83256) e 3 
+ 12(5 ^_ 4)7 (N + 4)(48C(3) (22200iV 6 + 11058CW 5 - 754767iV 4 

+ 13268217V 3 - 887361N 2 + 132711iV + 49508) + 144C(4)(5A r - 4) 3 

x (967V 3 - 276iV 2 + 221iV - 37) - 480C(5) (l5200iV 6 

- 79920iV 5 + 168488iV 4 - 183439iV 3 + 109827iV 2 - 34792iV 

+ 4656) - 99250iV 6 + 26925757V 5 - 36725295/V 4 + 140337844/V 3 

- 2306495707V 2 + 174742836/V - 50310868)e 4 , 
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1 



5JV-4 



(iV + 4)(iV-2)£ + 



(57V-4) 3 



(7V + 4)(30iV 3 + 257V 2 



- 217 N + 166) e 2 + 4{4 _\ N)5 (N + 4) (96C(3)(5iV - 4)(8iV 4 + 167V 3 

- 88iV 2 + 757V - 9) - 10307V 5 - 2751iV 4 - 46033iV 3 
+ 2(l03795iV 2 - 1336687V + 54904))e 3 

+ m l 5N)7 (N + 4) (48C(3) (l0300iV 7 - 39580iV 6 + 3284677V 5 

- 1208826/V 4 + 1806293iV 3 - 10742047V 2 + 94274/V + 82968) 
+ 144C(4) (2000iV 7 + 5200iV 6 - 58660iV 5 + 142451 TV 4 

- 160415iV 3 + 92612iV 2 - 25680/V + 2496) - 480C(5) (3800iV 7 

- 2280iV 6 - 47048iV 5 + 134947/V 4 - 1624217V 3 + 101278iV 2 

- 32704/V + 4448) - 39750iV 7 + 1242425iV 6 - 3090975/V 5 

- 35240910iV 4 + 171686590/V 3 - 295865304iV 2 
+ 226373292 TV - 65077096V 4 ■ 



The stability properties of the FPs are controlled by the eigenvalue exponents uoi of the 
stability matrix 

C df3 u (u,v) df3 u (u,v) \ 
d/3%,v) df3v("u,v) J (13) 
du dv ' 

taken at the given FPs. We calculated the eigenvalue exponents at all the FPs. For the most 
intriguing Bose and complex cubic FPs they are |3( 
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—e' + - 
20 



lr 257 384 



C(3) 



+ 



1 

2 £ ' 20" ' 8L 125 125 
Lr5109 + 624 ((3) _ 576 f(4) + 36« c(5) 



16 



1250 125 



125 



125 



co 2 



at the Bose FP and 
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To £ ~ Too 



1 



311 768 



C(3) 



1 

+ 16 



625 625 
« + &) + ^(4)-^<(5) 



250 3125 



625 



625 



1 

--£ 

2 



(60iV 3 - 1607V 2 + 1817V - 85 
4(57V-4) 2 (27V- 1) 



16(27V - 1)3( 5 JV - 4)4 l( 30720N7 ~ 178176iV6 + 456960iV5 
6483847V 4 + 5343367V 3 - 2519527V 2 + 62640 TV - 6336) C(3) 



(14) 
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+ 205607V 7 - 1653287V 6 + 644392iV 5 - 14068647V 4 

+ 17567457V 3 - 1224341iV 2 + 4337047V - 59052] e 3 

+ — 1- - [(99840007V 11 - 1492377607V 10 

64(27V - 1) 5 (57V - 4) 6 Lv 

+ 9211898887V 9 - 3096268032iV 8 + 63628458247V 7 

- 8473037952iV 6 + 7511188512iV 5 - 44527285927V 4 
+ 1734564864iV 3 - 4232962087V 2 + 581875207V 

- 3401472) C(3) + (-9216000 iV 11 + 77414400iV 10 

- 2990131207V 9 + 693626112iV 8 - 10644549127V 7 
+ 11278344967V 6 - 838693440 TV 5 + 4368034087V 4 

- 1559913127V 3 + 363703687V 2 - 4983552iV 

+ 304128)C(4) + (583680007V 11 - 5019648007V 10 

+ 1976709120iV 9 - 46536038407V 8 + 72243148807V 7 

- 77352729607V 6 + 5820821760iV 5 - 30785750407V 4 
+ 1122437280iV 3 - 2689905607V 2 + 381811207V 

- 2434560) C(5) + 81744007V 11 - 1116387207V 10 

+ 7829996007V 9 - 3263360352iV 8 + 85653920767V 7 

- 147036261727V 6 + 16826740917iV 5 - 128397201127V 4 
+ 63963919327V 3 - 1976965857iV 2 + 3404109167V 



24804436 



£ 4 



(TV — 2) (TV -l)(l407V 3 - 6007V 2 + 5677V -70) 2 
~2(57V-4) £ + 4(27V- l)(57V-4) 3 6 

- ^N-mL-A)* K 61440iv7 " 402432Ar6 + 975744iy5 

1 193856iV 4 + 8278087V 3 - 3375367V 2 + 77424iV - 7872) C(3) 

- 24880iV 7 + 1405447V 6 - 5509207V 5 + 15549287V 4 

- 25001557V 3 + 20908507V 2 - 8286287V + 1 18872] e 3 

~ WTTWTr ^^/rL JvF [(120064007V 11 - 892518407V 1 

64(27V- 1) 5 (57V -4) 7 Lv 

+ 4761141767V 9 - 18217536007V 8 + 43434487687V 7 

- 63730763527V 6 + 58385382247V 5 - 33253871527V 4 
+ 11263291687V 3 - 1980207687V 2 + 100468167V 
+ 986752) C(3) + (184320007V 11 - 1686528007V 10 
+ 6525158407V 9 - 14390231047V 8 + 20369387527V 7 



tW 
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- 19625103367V 6 + 13241966407V 5 - 6304728967V 4 
+ 209021904iV 3 - 462384007V 2 + 61724167V 

- 377856) C(4) + (-77824000iV 11 + 6790144007V 10 

- 24745369607V 9 + 50432870407V 8 - 64050662407V 7 
+ 52775692807V 6 - 28005888007V 5 + 8837846407V 4 

- 1 13276640iV 3 - 208950407V 2 + 95180807V 

- 1008640)C(5) - 24400007V 11 - 836416007V 10 

+ 10643847207V 9 - 55092820327V 8 + 162083669167V 7 

- 299180873527V 6 + 359181408997V 5 - 281888175157V 4 
+ 141727846897V 3 - 43331859347V 2 + 7212229167V 



49452712 



(15) 



at the complex cubic one. 

It is known that RG series are at best asymptotic, the coefficients of the series J2k fk9 h 
behave at large k as k\k b (—a) k [HJ, |32|. An appropriate resummation procedure has 
to be applied, in order to extract reliable physical information from them. To obtain the 
eigenvalues estimates, we have used an approach based on the Borel-Leroy transformation 

F(s; a, b) = ± M X) f e -ft Q* d Q W 



modified with the conformal mapping z = /^_\ + \ [25 , 26|, which does not require the 



knowledge of the exact asymptotic high-order behavior of the series [28]. The coefficients 



Ak(\) are determined from the equality B(x(z)) = where the Borel-Leroy transform 

B(x) is the analytical continuation of the series J2k a k v{l+k+i) xk absolutely convergent in the 
unit circle, are the coefficients of the original series. The numbers a and b characterize 
the main divergent part of the series. Since, in practice, we deal with a piece of the series 
only where the asymptotic regime might not be established, we vary parameters a and b in a 
neighborhood of their exact values. Our principle observation is that the result of processing 
F(e; a, b) exhibits very weak dependence on the transformation parameters a and b varying in 
a wide range. Thus, we put the stability of the result of processing with respect to variation 
of a and b into the foundation of our approach to resummation of divergent series^ [28 



The parameter A is chosen from the condition of the most rapid convergence of the series 
(|i~6[ ), that is from minimizing the quantity [1 — ^ | , where I is the step of truncation 
and Fi(e; a, b) is the /-partial sum for F{s\ a, b). If both eigenvalues are negative, the FP is 
infrared stable and the critical behavior of experimental systems undergoing second-order 
phase transitions is determined only by that stable point. For the Bose and the complex 
cubic FPs, our numerical results are displayed in Table I. 



3 The exact values of the asymptotic parameters a and b were calculated for the 0(n)-symmetric models 
, for the real cubic model pTJ , and recently for a number of more complicated models [R6( . Calculations 
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Table I: Eigenvalue estimates obtained for the Bose (BFP) and the complex cubic (CCFP) 
fixed points at N = 2 and iV = 3 in the four- loop order in e (e = 1) using the Borel-Leroy 
transformation with a conformal mapping. 



Type 
of FP 


N = 2 


N = 3 




UJ 2 




UJ 2 


BFP 


-0.395(25) 


0.004(5) 


-0.395(25) 


0.004(5) 


CCFP 


-0.392(30) 


-0.029(20) 


-0.400(30) 


-0.015(6) 



It is seen that the complex cubic FP is absolutely stable in D — 3 (e = 1), while the Bose 
point appears to be of the "saddle" typef]. However ^'s of either points are very small at the 
four-loop level, thus implying that these points may swap their stability in the next order 
of RG approximation. We can compare u 2 at the complex cubic FP quoted in Table I with 
the three-loop results of Ref. [ID] obtained in the framework of RG approach directly in 3D. 
Those estimates uo 2 = —0.010 for N = 2 and u 2 = —0.011 for N = 3 are solidly consistent 
with ours. 

In addition to the eigenvalues, we have computed the critical dimensionality of the order 
parameter of model (§). The critical spin dimensionality iV c for the complex cubic model 
is determined as that value of N, at which the complex cubic fixed point coincides with 
the isotropic one. Equivalently, for N = N c the second eigenvalue of the stability matrix Q 
taking at the complex cubic FP vanishes, uj 2 = 0. The four-loop e expansion reads 



Nr. 



2 " £+ 24 



6C(3) - 1 



s 2 + 



144 



45C(3) + 135C(4) - 600C(5) - 1 



Instead of processing this expression numerically, we have established the exact relation 
N c = 2 n c between the critical spin dimensionalities of the complex and the real hypercubic 
models, which is independent on the order of approximation used. In fact, both models, 
the complex and the real cubic ones, exhibit effectively the isotropic critical behavior at 
N = N c and n = n c , respectively. Therefore, because the complex 0(2iV)-symmetric model 
is equivalent to the real 0(n)-symmetric one, the relation 2N C = n c holds true. For N > N c 
the complex cubic FP of model (|2|) should be stable in 3D. 

The five-loop e expansion for n c was recently obtained in Ref. |33]]. Resummation of 
that series gave the estimate n c = 2.894(40) (see Ref. [ff^fl). Therefore we conclude that 
N c = 1.447(20) from the five-loops. Practically the same estimate N c = 1.435(25) follows 
from a constrained analysis of n c taking into account n c = 2 in two dimensions P7[ . From 



the recent pseudo-e expansion analysis of the real hypercubic model [55| one can extract 
N c = 1.431(3). However the most accurate estimate N c = 1.445(20) results from the value 
n c = 2.89(4) obtained on the basis of the numerical analysis of the four-loop |34] and the 



show that the resummation method employed in this paper gives the same results as the conventional 
technique using the exact values of the asymptotic parameters (see Rcfs. 27 and p8[) 



9 The fixed point is called to be of the "saddle" type provided their eigenvalue exponents oji and U2 are 
of opposite sings at the (u, v) plane. 
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six-loop [^7| 3D RG expansions for the /3-functions of the real hypercubic model. Since 
N c < 2, the critical thermodynamics in the NbO"2 crystal, in the antiferromagnets TbAu2, 
DyC2, K 2 IrCl6, TbD 2 , MnS2, and Nd as well as in the gelical magnets Tb, Dy, and Ho should 
be controlled by the complex cubic fixed point with a specific set of critical exponents, in 
the frame of the given approximation. 



Ill Critical Exponents and Conclusion 

In the previous section we have shown that one of the four FPs is absolute stable in 3D. 
This is the anisotropic complex cubic FP. Now let us calculate the critical exponents. To do 
this, substitute the coordinates of the FPs found in Sec. II into the expressions for 77 and 
7 _1 [see Eqs. (|10| ) and (|12j)]. For the stable complex cubic FP, this yields 



4(57V-4) 2 
5357V 2 + 6527V 



(TV — 1)(2N — 1) + 



16(57V-4) 4 



(N-l) 



19(W d 



324 



64(5iV - 4) 6 
768iV 4 + 15727V 3 - 1613iV 2 + 777 iV - 132 
226917V 4 + 685277V 3 - 2803997V 2 + 3268887V - 127676 



96C(3)(l607V 5 
105707V 5 



(17) 



7" 



3 e e 2 
2(57V-4)^ ~ ' + 4(57V-4) 3 



(N-l) 



407V 2 - 2147V 



205 



8(57V-4) 5 
1597V- 13 



(TV 



12C(3)(57V-4)(327V 3 - 1567V 2 
671027V 



32558 



9407V 4 - 67487V 3 + 426817V' 

(TV - 1) 8C(3) (57V - 4) (56507V 5 



65510 



64(5TV-4) 7 

+ 58655TV 4 - 303683TV 3 + 396967TV 2 - 965027V - 
+ 72C(4)(5TV - 4) 3 (32TV 3 - 156TV 2 + 1597V - 13) 
- 160C(5) (57V - 4) 2 (3047V 4 - 16167V 3 + 24247V 2 - 12777V 
+ 265) - 543007V 6 - 172207V 5 - 74228277V 4 + 424275647V 3 



- 865211377V 2 + 765634487V - 25005620 

The other critical exponents can be found through the known scaling relations. The nu- 
merical estimates obtained via transformation ( |T5] ) are collected in Table II. The critical 
exponents for the isotropic and Bose FPs are also presented, for comparison. We can com- 
pare our results with the available experimental data. For example, in the case of the 
structural transition in Nb0 2 crystal, the critical exponent of spontaneous polarization was 
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Table II: Critical exponents for the isotropic (IFP), the Bose (BFP), and the complex cubic 
(CCFP) FPs at N = 2 and N = 3 calculated in the four-loop order in e (e = 1) using the 
Borel-Leroy transformation with a conformal mapping. 



Type 
of FP 


N = 2 


N — 3 


7] 


V 


7 


V 


V 


7 


IFP 


0.0343(20) 


0.725(15) 


1.429(20) 


0.0317(10) 


0.775(15) 


1.524(25) 


BFP 


0.0348(10) 


0.664(7) 


1.309(10) 


0.0348(10) 


0.664(7) 


1.309(10) 


CCFP 


0.0343(20) 


0.715(10) 


1.404(25) 


0.0345(15) 


0.702(10) 


1.390(25) 



measured in Ref. ||, 0.33 < (3 < 0.44. Our estimate (3 = 0.371 obtained using the data of 
Table II and scaling relations lies in that interval. 

In conclusion, the four-loop e-expansion analysis of the Ginzburg-Landau model with 
cubic anisotropy and complex vector order parameter relevant to the phase transitions in 
certain antiferromagnets with complicated ordering has been carried out. Investigation of 
the global structure of RG flows for the physically significant cases N = 2 and N = 3 
has shown that the anisotropic complex cubic FP is absolutely stable in 3D. Therefore the 
critical thermodynamics of the phase transitions of concern is governed by this stable point 
with specific critical exponents. The critical dimensionality N c = 1.445(20) obtained from 
six loops supports this conclusion. At the complex cubic FP, the critical exponents were 
calculated using the Borel-Leroy resummation technique in combination with a conformal 
mapping. For the structural phase transition in Nb02, for the antiferromagnetic phase 
transitions in TbAu2 and DyC2 as well as for the phase transitions in the rare-earth metals 
Ho, Dy, and Tb, they were shown to be close to the critical exponents of the 0(4)-symmetric 
model. In contrast to this, the critical exponents for the antiferromagnetic phase transitions 
in K 2 IrCl 6 , TbD 2 , MnS 2 , and Nd turned out to be close to the Bose ones. 

Although our calculations show that the complex cubic FP, rather than the Bose one, is 
stable at the four-loop level, the eigenvalues \u 2 \ of both points are very small. Therefore 
the situation is very close to marginal, and the FPs might change their stability to opposite 
in the next order of perturbation theory, so that the Bose point would occur stable. This 
conjecture is in agreement with the recent six-loop RG study of three-coupling-constant 
model ([!]) directly in three dimensions [ 3~6"| . There, the authors argue the global stability 



of the Bose FP, although the numerical estimate u 2 = —0.007(8) of the smallest stability 
matrix eigenvalue at the Bose FP appears to be very small and the apparent accuracy of 
the analysis does not exclude the opposite sign for u 2 - In this situation, it would be highly 
desirable to compare the critical exponents values obtained theoretically with values that 
could be determined from experiments, in order to verdict which of the two FPs is really 
stable in physical space. Finally, it would be also useful to investigate certain universal 
amplitude ratios of the model because they vary much more among different universality 
classes than exponents do and might be more effective as a diagnostic tool. 
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